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A THEORETICAL INVESTIGATION OF THE INPUT CHARACTERISTICS 


OF A RECTANGULAR CAVITY-BACKED SLOT ANTENNA* 

C. R. Cockrell 
Langley Research Center 

SUMMARY 

Equations which represent the magnetic and electric stored energies are derived 
for an infinite section of rectangular waveguide and a rectangular cavity. These repre- 
sentations which are referred to as being "physically observable" are obtained by con- 
sidering the difference in the volume integrals appearing in the complex Poynting theorem. 
It is shown that»the "physically observable" stored energies are determined by the field 
components that vanish in a reference plane outside the aperture. 

These "physically observable" representations are used to compute the input admit- 
tance of a rectangular cavity-backed slot antenna in which a single propagating wave is 
assumed to exist in the cavity. The slot is excited by a voltage source connected across 
its center; a sinusoidal distribution is assumed in the slot. Input-admittance calculations 
are compared with measured data. In addition, input-admittance curves as a function of 
electrical slot length are presented for several size cavities. 

For the rectangular cavity-backed slot antenna, the quality factor and relative band- 
width were computed independently by using these energy relationships. It is shown that 
the asymptotic relationship which is usually assumed to exist between the quality band- 
with and the reciprocal of relative bandwidth is equally valid for the rectangular cavity- 
backed slot antenna. 


INTRODUCTION 

The aperture (or slot) antenna is one of the most widely used antennas because it is 
relatively simple to build and can be flush mounted in conducting bodies such as in the 
surface of aircraft and spacecraft, thus becoming an integral part of the vehicle. For 
such applications the aperture antenna also meets the requirements of small size and low 
weight. In addition to its practical usefulness, radiation and impedance characteristics 
can be investigated theoretically without too much difficulty (refs. 1 to 3). Such inveSti- 

*The information presented herein was offered as a dissertation in partial fulfill- 
ment of the requirements for the Degree of Doctor of Philosophy in Electrical Engineer- 
ing, North Carolina State University, Raleigh, North Carolina, May 1974. 


gations are usually conducted first for an ideal model such as a narrow slot in a perfect 
conductor of infinite extent. 

The input admittance Y s (or impedance) of a narrow slot in a perfectly conducting 
infinite sheet can be determined via the Booker relationship (ref. 4) whenever the slot is 
free to radiate on both sides of the infinite sheet. This relationship, which can be found 
in many books (refs. 5 to 8), is given by Y s = AZ^jz^ where Z^ is the input imped- 
ance of the complementary dipole (planar dipole) and Z Q is the characteristic imped- 
ance of the surrounding medium. In practical applications the slot is backed by some 
sort of cavity, thus destroying the symmetry upon which the Booker relationship depends. 
The cavity-backed aperture antenna has been the subject of many papers over the past 
two decades, or longer (refs. 9 to 13). 

When the slot is backed by a cavity on one side of the infinite sheet, the radiation 
pattern and impedance characteristics of the slot antenna are altered; the radiation 
resistance, the bandwidth, and the stored energy are changed. The impedance properties 
of the apertures (and slots) which are backed by a rectangular cavity have been investi- 
gated by many authors (refs. 9 to 13). In references 9 and 10 the backing rectangular 
cavity was a shorted waveguide whose cross section was the same as that of the aperture; 
reference 9 is further restricted to small cavities. In references 11 to 13 the thin slot is 
backed by a rectangular cavity of different cross section. 

In references 9 and 10 the relationship between quality factor and the reciprocal of 
relative bandwidth, known to exist in nonradiating systems, is assumed equally valid for 
cavity-backed aperture antennas. Quite often this relationship is assumed for antennas 
in general (refs. 14 and 15). In reference 16, quality factor and inverse bandwidth are 
related in an order of magnitude sense. The antenna in references 14 and 15 is assumed 
to be such that its aperture distribution is frequency independent; hence the frequency 
derivative of its reactance or susceptance is shown to be proportional to the total stored 
energies. The reciprocal of relative bandwidth would then be equal to quality factor. 

For planar antennas in which the aperture distributions are frequency dependent, the 
reciprocal of relative bandwidth is proportional to the frequency derivative of the differ- 
ence in stored electric and magnetic energies; whereas the quality factor is proportional 
to the total stored energy (ref. 17). The "redefined” definition of quality factor given by 
Rhodes in reference 18 should be used. 

The establishment of a relationship between the quality factor and the reciprocal of 
relative bandwidth for the aperture antenna would be of analytical importance in the area 
of antenna synthesis (refs. 17 and 18). The evaluation of quality factor is at a single fre- 
quency (resonant frequency); whereas the evaluation of relative bandwidth requires a 
knowledge of the complete frequency behavior of the antenna. The quality factor and 
relative bandwidth for the planar dipole antenna have been shown to be related in a recip- 
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rocal manner by Rhodes (ref. 17). He established this relationship by calculating the 
quality factor and relative bandwidth from independent equations which were based on 
his time-average "physically observable" stored-energy representations. By consider- 
ing the difference of the volume integrals that appear in the complex Poynting theorem, 
Rhodes was able to show that the infinities associated with the individual volume integrals 
canceled exactly, leaving what he refers to as time-average "physically observable" 
stored magnetic and electric energies. 

In the present paper, expressions which represent the input admittance, the quality 
factor, and the relative bandwidth of the rectangular cavity-backed slot antenna are 
derived. For an assumed sinusoidal slot distribution and a single propagating wave in 
the cavity, input-admittance calculations are compared with available measured data. In 
addition, input-admittance curves as a function of electrical slot length are given for 
several size cavities. The primary purpose of this paper, however, is to determine 
whether the method of Rhodes implies a reciprocal relationship between the quality fac- 
tor and relative bandwidth for a rectangular cavity-backed slot antenna. Thus Rhodes' 
concept of time-average "physically observable" stored energies is used in determining 
this relationship. A pair of time-average "physically observable" stored energies for the 
internal part of the antenna (cavity or infinite waveguide) are derived in terms of ampli- 
tude coefficients at a reference plane. These coefficients are related to the assumed 
sinusoidal distribution in the narrow slot aperture by applying the appropriate boundary 
conditions. The cross sectional dimensions of the cavity are chosen so that only one 
propagating wave exists. The depth of the cavity is assumed to be deep enough so that 
its back wall will not influence the assumed slot distribution. 

The external part of the antenna is the half-space. The solution to this part is 
given by Rhodes in reference 17 for a slot in an infinite ground plane. 

The internal and external parts of the solution are then combined by applying the 
complex Poynting theorem to the small volume which is formed by the thickness and the 
openings of the slot. This volume encloses a voltage source which is applied across the 
center of the slot. As the slot thickness shrinks to zero, the input power equals the sum 
of the power which enters the cavity (finite volume) and that which enters the half-space 
(infinite volume). Once the internal and external parts have been combined, the input 
admittance, the quality factor, and the relative bandwidth are computed numerically for a 
number of cavity and slot sizes. The only power loss is from radiation since all the 
metallic surfaces are assumed to be perfect conductors and the region inside the cavity 
is assumed to be a vacuum. 


3 



SYMBOLS 


A mn > B mn> 

l 

Cmn)Dmn^ 

> amplitude coefficients 

Ai,A2 

slot aperture (see fig. 2) 

a 

height of waveguide (or cavity) 

a* 

width of slot (see fig. 1) 

B 

input susceptance 

B.W. 

relative bandwidth 

b 

width of waveguide (or cavity) 

Cin(x) 

. (* 1 - cos u du 
J 0 u 

d 

depth of waveguide (or cavity) 

E 

electric field 

E xa» E ya 

aperture electric field 


f(y) Fourier transform of S(u>) 

G input conductance 

H magnetic field 

I current 


J 

j 

K()(ay) 


current density 

modified Bessel function of the second kind 


4 



k 

kx’Vkz 

a 

£oA 

m,n 

p i 

Q 

Si(x) 

S(w) 

s l, s 2 

t 

V 

v 0 

v 

v' 


«W e » 


«w m » 


x,y,z 


= 

directional wave numbers 
length of slot (see fig. 1) 
resonant slot length 
modal numbers 
power loss in cavity 
radiated power 
quality factor 

nX . 

= f 

J 0 U 

Poisson sum 

surface 

time, sec 

voltage 

voltage across center of slot 
volume of cavity 

volume of semi-infinite free space 
volume of slot 

"physically observable" electric stored energy 
"physically observable" magnetic stored energy 
Cartesian coordinates of rectangular waveguide 
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x’,y',z' Cartesian coordinates of slot 

Y input admittance 

Y g input admittance of slot antenna 

Z d input impedance of complementary dipole antenna 

Z Q characteristic impedance of free space 

6 depth of slot (see fig. 1) 

e permittivity of free space 

X wavelength 

/u. permeability of free space 

(*) angular frequency 

oo T angular resonant frequency 

Subscripts: 

ap aperture 

inf infinite waveguide 

non nonpropagating wave 

pro propagating wave 

Asterisks denote complex conjugates. 

GENERAL EXPRESSIONS RELATING INPUT ADMITTANCE, 

QUALITY FACTOR, AND RELATIVE BANDWIDTH 
TO STORED ENERGIES 

The rectangular cavity-backed slot antenna is shown in figure 1. The antenna is 
divided into two parts: internal and external. The external part of the solution has 
already been solved in reference 17. The internal part of the solution is solved in the 
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X 



x f 



Figure 1.- Slot backed by a Figure 2.- Slot geometry. (Aperture A 1 

rectangular waveguide. empties into half-space and aperture Ag 

empties into the rectangular cavity. ) 

present paper by using the concept of time-average "physically observable" magnetic and 
electric stored energies. 

However, before determining these energies, the internal and external parts of the 
solution are combined by applying the complex Poynting theorem to the slot volume shown 
in figure 2. (See ref. 5.) Thus, 

10 Ex5 *'dS=-i2“feOTS.H*dv-|yj'j'E.E*dv) -ij'j'j'E.r dv (1) 

V v" v" / v M 


4OT 


J* dv 


= | EXH* 
2 Ai 


(-z dx’ dy’) +| 

a 2 


E x H* • (+z dx' dy’) 




( 2 ) 


The term 


4OT 


E • J* dv 


represents input power; therefore, 
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( 3 ) 


" 5 Iff E • J* dv = +| W* 

V” 

The first two terms on the right-hand side of equation (2) are the complex power flow into 
the half-space and cavity, respectively, and the last term represents the net time-average 
stored energy in volume v". 

As 6 approaches zero, the volume v" approaches zero, and, hence, the net 
stored energy in v" is zero. Therefore, 

| VI * = Ijj' E x H* • (-z dx’ dy') + 1 E x H* . (+z dx’ dy’) (4) 

Ai A 2 

The complex Poynting theorem will now be applied to the volume v* which is infi- 
nite and to the volume v which is finite (cavity). (See fig. 1.) The complex power flow 
into the half-space z < 0 is given by 


E XH* • (-z dx' dy’) = P r + j2a jJJ H • H* dv - | JJJ E ■ E* 
Af \ v’ v’ 



( 5 ) 


and the complex power flow into the cavity z > 0 is given by 


JJexh*. 

■ (+z dx' dy') = P z + j2o> JJJ H • H* dv 


Ai 

V v 

v / 


( 6 ) 


where P r is the radiated power and P^ is the power loss in the cavity. The volume 
integral terms represent the time-average net stored energies in their respective 
volumes. . 

Assuming the power loss P^ in the cavity is zero, equations (4), (5), and (6) can 
be written as 
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Now, let I = YV, so that 



Taking the conjugate, equation (8) becomes 



The admittance is now given as 
Y - G + jB 



— v*v 

2 


( 7 ) 


( 8 ) 


(9) 


( 10 ) 
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Therefore, 


G = 


- W* 


(ID 


B = 


2 “(f III 5 ■ S * dT - fill E ' r dV ) 2l 4lII 3 • S* dv - E • 

\ V* V* j \ V V 


E* dv] 


i vv * 




B 


2w^«W m » v t - «W e » v ') 2oi^«W m » v - «W e » v 


1 VV* 


I w* 


(12a) 


(12b) 


B = 


-2a>(«W m » - «W e ») 


J: VV* 

2 


(12c) 


where 


<< W m >> - «W e » = ^«W m »v' - «W e »v') + («W 


m » v - «W e»v 


For a high Q system with a resonance frequency of u> r , 


4G -0 


dw 


(i), 


B(co r ) = 0 


As one moves off resonance, the input admittance can be written as 


(13) 


Y(o>) *G(u> r ) + j(u> - w r )^j 


(14) 
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At the half-power point, 


( \ j — j “ 

co r + = |/2 Y(co r )e 4 = ]/2 G(co r )e 4 


( 15 ) 


Substituting equation (15) into equation (14) and equating imaginary parts gives 


(co - co r )| 


(— ) 

\ dco /co=CO T 


~ G(co r ) 


(16a) 


or 




(16b) 


Therefore, the relative bandwidth is 


B.W. = 


A co 2G(co r ) 


dB 

dco 


co T 


CO r 


(17) 


Multiplying both numerator and denominator by ^ VV* gives 


B.W. = 


i VV* G<o> r) ] 


2P. 


co r I VV* — co r — VV* — 

r 2 dco r 2 dco 


(18) 


From equation (12c), 


dB _ _d_ 
dco dco 


-2co(«W m » - «W e ») 


— VV* 
2 


so that 


1 d 


dB\ 

da Va>=co r A VV* 
2 


-2co(«W m » - «W e ») 


co=co r 


(19) 
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( 20 ) 



z, \ r 

t \ 

I**) = .i 2 

Ww= o, r dwl - 

wK < W m » - «W e »J 


U>=U) r 


and, thus. 


B.W. = 


2P t 


-co r 2a>(«W m » - «W e ») 




(21a) 


or 


B.W. = 


2P t 


-U) 



2o)^«W m » v t - «W e » v ') + -p 2w^«W m »v - «W e » v j 


(21b) 


0>=Ct) T 


The concept of quality factor Q is usually defined as being proportional to the sum 
of the time-average magnetic and electric stored energies. However, an inconsistency in 
this classical definition has been asserted by Rhodes in reference 18. This definition 
includes parts of the magnetic and electric stored energies that can not be observed since 
these energies are obtained by treating the volume integrals individually in the complex 
Poynting theorem. This inconsistency has been removed by Rhodes by redefining Q in 
terms of the time-average ’’physically observable" stored magnetic and electric energies, 
determined through differencing the volume integrals in the complex Poynting theorem. 
This "redefined" definition of Q is 


Q = 


0)( << W m >> + «W e ») 


UJ = U)y 


( 22 ) 


In terms of the volumes considered in the present paper, the quality factor becomes 


Q = 




(jO — LOy 


(23) 


12 


J 



But at resonance the observable magnetic and electric energies are equal, that is, 


«W jjj T + «W m »v — «W 0^y* + ^W e ^>v 


( 24 ) 


Hence, 


Q = 


2w(«W m>e » v » + << W m>e » v ^ 


u>=w T 


(25) 


where the subscript m,e means that either magnetic or electric "physically observable" 
stored energy is used. 


DERIVATION OF TIME-AVERAGE MAGNETIC AND ELECTRIC VOLUME 
INTEGRALS IN A RECTANGULAR WAVEGUIDE SECTION 


The volume integrals appearing in the complex Poynting theorem are given as 
tL H • H* dv and - E • E* dv. The fields for a rectangular waveguide will now 
v v 

be substituted into each of these volume integrals. Each volume integral will be carried 
through in its entirety. 

The total magnetic field is 


H = xH x + yHy + zH z 


(26) 


where the components are given by equations (A23) in appendix A. Substituting equa- 
tion (26) into the volume integral results in 


^ H • H* dv = §§§ ^xH x + yH y + zH z ) • (xH* + yH* + zH*)dv (27) 

v v 

where the volume is bounded by the walls of the rectangular waveguide and transverse 
planes Si(where z = 6) and S 2 (where z = d). Therefore, 

£ H • H* dv = i±- (h x H x + H y Hy + H z H z )dv (28) 

v v 
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Each term of the integrand (eq. (28)) is written out explicitly as follows: 



The integration over the cross section is performed first. Due to orthogonality, the inte- 
grations for m * m' and n * n' are zero. In addition, the resulting integrations are 
weighted differently depending on m and n. When m = 0 in equation (29a) and n = 0 
in equation (29b), the integrations are identically zero. When m * 0 and n = 0 in 
equation (29a), n * 0 and m = 0 in equation (29b), and m * 0, n = 0, n * 0, and 
m = 0 in equation (29c), the integrations result in a multiplying factor of 1/2. When 
m 4 - 0 and n * 0, the integrations of equations (29) produce a multiplying factor of 1/4. 
For m = 0 and n = 0, the multiplying factor resulting from the integration of equa- 
tion (29b) is 1. The integrations over the cross section for equations (29) can, therefore, 
be written as 
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X io "’" 5 " ^ = “ Jo - {[>? * (“/]- 

* (?)f^--}{[-5 *(#-■ - {£ * f?)j»~ • ‘W’K 

* W0 c ™}'’ (> " lJ ’ ' {E 5 * (“)’]•- • P® “"){}’ * (fi ! ]“'"- •Pj!?) t -'-}‘ , ’*^ ! ' 


3 I H z Hjdxdy=f l l iEfa 

y=0 x=0 ' “ 


(™1 Cmn - HZL A )(™c* -HI A* \ e ’ j ( kz ' k z) Z + (mv r nv A \/mv ott n , \-j(*z+kJ)2 

U mn b A mn^ a C mn A mn je + c mn - — A mn j^— D mn - — B mn je 

Dmn - — B ](™1 C* -Ha* \j( k z +k 3 z /imr n* R \(mv n * r»7 n * \ + j( k z‘ k £) z | 

\ a Umn b mn J^~ Cmn V m T + \^ mn "T mn J(lT mn “T mn J 


c m,n 


(m,n = 0) 
(m,n * 0) 

(m,n = 0) 
(m,n ± 0) 


Now, rewrite equations (30) as 


j I H X „>d»=f x Z^Zti 

y=0 x=0 m=0 n=0 ^z^z w M 


-i(kz-^)z 


{h 2 )(w t ) A — ♦ [(rf * ^ c -} ((” x? h 


+ 1 

(nii\ 2 ,2 


(nni\(nit\ n 

/nr\ 2 . 2 



(t) +kz 

D miy e j + 


(“j ^ 

D mn 


j(^z _k z) 3 


J I H y H J dx dy = 2^ £ £ 

y=0 x=0 m=0 n=0 f 1 





r\ 2 l /r 


-jfkz-k^z // 

♦(? 

■j A mn + (- 

r)(f) c ™ 

C 

- 2R \\ 





2 

z\ 

2 /mir\ 2 

j \ 

/ \ 

j(k z -k3z 

) + 



[f) D ™ 

e 


(f)] A ”*(f)(T) c "){[^(ff 


(30b) 


(30c) 


(31a) 


(31b) 


(32a) 


(32b) 
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(32c) 


b a 


y J H z H’dxdy=^ 

y=0 x=0 




~j(^~kz) z 

e 


+ 2 Re 


{?< 




Before integrating on z, the wave number kg will be written explicitly for propagating 
waves as 



and for nonpropagating waves as 




(33a) 



Substituting for k z and integrating on z (the integration on z is over the limits 
z = 6 and z = d; the distance 6 is then allowed to approach zero), the total magnetic 
volume integral becomes the following: For propagating waves, 



(Equation continued on next page) 
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d + 2 Re 




( 34 ) 


For nonpropagating waves, 



(Equation continued on next page) 
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The total magnetic volume integral is the sum of the volume integrals for propagating and 
nonpropagating waves; that is, 


OT fia *H£ftI fl - rdv 




H • H dv 


Attention is now turned to the volume integral of the electric field. The total elec- 
tric field is given as 

E = xE x + yEy + zEz (37 

where the components E x , Ey, and E z are given by equations (A18). Substituting 
equation (37) into the volume integral results in 




4 J J J (E X Ex + EyEy + E Z E Z dv 


i^x = £ ^ (A mn e jkzZ + B mn e jkzZ ) cos ~ x sin ^ y 

m=0 n=0 


00 OO , 

1 1 ( A m'n' 


m'=0 n'=0 


Jk z z * ]k z z\ m' 7 r . n' 7 t 

e + B . „,e cos x sin 

m'n' / a b 
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E y E y - 


00 oo 

1 I 

m=0 n=0 


rmn e 


" j ^ Z + D mn e jkz jsinB£xcos^y 


OO oo 

2 I 

m=0 n=0 


(*-'m'n ,e 


+jk*z 


z 

+ D m'n ,e 


-K z 


z . m ? 7 r n f 7r 

/sin x cos y 

' a b 


(39b) 


oo oo 


"■"HI 


m=0 n=0 


-jk z z 


mir 


B 


mn 


f D mn) 


jkz z 


sin x sin ^ y' 
a b ' 


[ OO oo 

y y +— 

/mil a* 

n' 7 r r * ' 

\ jkz z / m t, 

|m’= 0 nM) J kz 

y a rn'n' 

b m'n; 

I 6 - 1 “ 


R * 

D TV.'r>' 


n’ir n * V 

y D m'n" 


\ "jk^z 

m T 7r . n T 7r 1 

r 

sin — — x sin - — y' 
a b J ( 


j 


(39c) 


As before, integration over the cross section is performed first, i.e., 

J | E-E*dxdy=^ £ J fvm( A mne' JkzZ +B mn e )kzZ )(A^ ne ,kJz 
y=0 x=0 m=0 n=0 L 


+ Bmn e 


•i * 

-jk z z 


j k z z R J k z z )L* Jk z z p.* -jk z z 


+ Vn| C mn e + D mn e j( C mne + D mn e 
. 'Wnl 


k z k z 
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where 



and orthogonality has been taken into account. 

The volume integrals of the electric field for propagating and nonpropagating waves 
will now be determined in a manner similar to that for the magnetic volume integrals. 
These electric volume integrals become the following: For propagating waves (eq. (33a)), 



For nonpropagating waves (eq. (33b)), 
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Therefore, the total electric volume integral is 



Derivation of Time-Average "Physically Observable" Stored Energies 
of an Infinite Rectangular Waveguide 

When the complex Poynting theorem is applied to radiating systems, or to any 

closed volume, the volume integrals of H • H* and i. E • E* appear in the equation 

4 4 

only through their difference. However, in the literature (e.g., refs. 5, 6, and 8) these 
volume integrals are usually interpreted individually. It has been shown by Rhodes in 
reference 17 and Collin and Rothschild in reference 14 that by considering an infinite 
volume (planar antenna) for certain aperture distributions such an interpretation leads to 
infinite stored energies. Rhodes showed that by treating the difference of the volume 
integrals the infinities cancelled exactly. From the remaining parts of the difference he 
determined integral representations of the "physically observable" time-average magnetic 
and electric stored energies. 

An analogous problem to the planar antenna would be the infinite rectangular wave- 
guide which can be identified with figure 1 for the case where d is infinite, a' = a, and 
£ = b. Applying the complex Poynting theorem to the volume v shown in figure 1 pro- 
duces the volume integrals given by equations (36) and (44). The surface S 2 at z = d 
is treated as a mathematical boundary; that is, the coefficients B m n and Dmn corre- 
sponding to no wave traveling in the zrdirection are zero. For a source (aperture) that 
is located at the z = 0 plane, with perfectly conducting walls and no losses within the 
waveguide, radiation occurs through z = d as d — °°. Thus, setting B mn and D mn 
equal to zero in equations (35) and (43) and allowing d to approach infinity, the magnetic 
and electric volume integrals become 
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J177 

+ 

b 




(46) 


It is shown in appendix B that the contributions to the electric and magnetic volume 
integrals from the propagating waves are identical. Therefore, when the volume integrals 
are differenced, these contributions cancel identically. In addition, it is also shown in 
appendix B that further cancellations occur in the volume integrals when the terms are 
written in a certain form. Before discussing these cancellations, it is essential to offer 
some physical interpretation which would justify writing the terms in a proper form for 
carrying out these cancellations. 

The coefficients A mn and C mn are determined by matching the aperture distri- 
bution with the waveguide fields at z = 6 as 6 tends to zero; that is, 


E xa (x,y,k) = lim E x (x,y,6,k) = lim Y Y A mn cos 211 x sin 22 y e * z (47a) 

6—0 6-0 Jr'n _ a a b 


m=0 n=0 


E ya (x,y,k) = lim E y (x,y,6,k) = lim ) Y C mn sin 222 x cos H y e ^ z6 (47b) 

6-0 6-0 a b 


m=0 n=0 


Multiply equation (47a) by cos 21! x sin 21 y and equation (47b) by sin 21Z x cos 21 y> 

3 D 3D 

and integrate over the waveguide cross section. For planar apertures in perfectly con- 
ducting metal, E xa and E ya outside the aperture are zero; hence, 

f f E xa cos 211 x sin 21 y dx dy 

Amn = — — 7 T (m,n * 0) (48a) 
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(48b) 


(48c) 


(48d) 


These coefficients are finite for all values of frequencies since the electric-field distri- 
butions in the aperture are well-defined functions of frequency. Taking this into account, 
the volume integrals given by equations (45) and (46) become infinite, in general, as 
k z — 0 because the contributions from the E z , H x , and Hy components to the volume 
integrals become infinite. The contributions from the E x , Ey, and H z components 
are finite. When the volume integrals are differenced, these infinities must cancel in 
some way since the surface integral in the complex Poynting theorem is finite for all fre- 
quencies. Following Rhodes (ref. 17), one may argue that since these canceling terms 
disappear from the complex Poynting theorem for all frequencies they must then have no 
physical significance and may therefore be neglected in the volume integrals, even when 
these integrals are considered separately. Those terms which do not so cancel are 
accordingly referred to as "physically observable." The difference is written as 
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It is shown in appendix C that the expressions associated with the two braced terms in 
equation (49) vanish. Therefore, the "physically observable" stored energies for an 
infinite rectangular waveguide are 


«W m » in f = 


ab 


ou 

I 1 


16w M m =0 n=0 


e m e n 


m7r 


'mn 


-S2LA 


mn 


llffif + Iff - k 2 
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k2 <(T)\f 2 


(50a) 


OO OO 




m=0 n=0 




nm 


l mn 


+ Vn 


'mn 


+ ^ -k< 

b 


k 2 jBlfjslf 




(50b) 


Notice the similarity in this result (eqs. (50)) to the result given by Rhodes for the planar 
antenna (ref. 17), which is another radiating system with infinite volume. The energies 
stored in a section of uniform waveguide are usually given in the literature as energy per 
unit length (see ref. 19). To determine the energies stored in a specified section of length 
one should then multiply by the length. However, for an infinite length such procedure 
would result in infinite stored energies. The stored energies represented by the pre- 
ceding equations (50) show that they remain finite even for an infinite section of waveguide. 


Derivation of Time- Average "Physically Observable" Stored Energies 

of a Rectangular Cavity 

When a perfect electric conductor is placed at z = d, the reflected amplitude coef- 
ficients B mn and D mn are not zero. Instead, these coefficients are related to the 
amplitude coefficients A mn and C mn , respectively. This relationship is established 
by applying the boundary condition at z = 0 and z = d on the transverse electric fields. 
At z = 0, the coefficients are determined in the same manner as given before; that is, 

§§ E xa cos x sin ^ y dS 

^mn + ®mn = ( m > n ^ 0) (51a) 

ab/4 
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^mn + ®mn - 


(m,n * 0) 


(51b) 


Jf EyasmiSIxcosHyds 

i±E ^ 

ab/4 

E xaSin2 IydS 

A 0n + B 0n = ( n * 0) (51c) 

ab/2 

^ Ey a sin ^IZL x dS 

C m 0 + B m0 = ^ ( m ^ 0) (5 Id) 

ab/2 


At z = d, 

-jk z d jk z d ^ 

A mn e + E mn e = 0 

> 

_ -i k zd jk z d 

Cmn e + B mn e - 0 

J 

Hence, 

“N 

-j2k z d 

B mn = "Amn e 

> 

D - -C e' J2kzd 
^mn “ L mn e 


(52) 


(53) 


where k z is given by equations (33). By substituting equations (53) for B mn and D mn 
into equations (34), (35), (42), and (43), the magnetic and electric volume integrals become 
the following: 
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By expressing the exponentials in their trigonometric and hyperbolic representations, 
equation (54) is written as follows: 



(Equation continued on next page) 
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When the volume integrals are differenced, cancellations in terms occur as in the infinite 
waveguide case. However, before discussing these cancellations, it is necessary to give 
some physical interpretation leading to these cancellations. 

Substituting equations (53) for Bnin 3-nd Dmn into equations (51) results in 


l mn - 


ft 

ap 


E xa cos x sin — y dx dy 


3^1 . e ' j2k z d 


ffap Eya sin ^ x cos f y dy 


'inn 


-4 - e- J2kad 


E xa sin 22 y (Lx dy 

Ao„>— ^ 


n - 


(58a) 


(58b) 


(58c) 
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^mO _ 


(58d) 



Ey a sin x dx dy 




For the case of the infinite waveguide, the application of Rhodes' method of identi- 
fying the "physically observable" stored energies has been seen to involve the exclusion 
of certain terms in the stored energy expressions which, when considered separately, 

tend toward infinity as l/k^ when k z — 0. For the finite-length waveguide (cavity) 

/ o ° 

such l/k z terms also exist, but they are not the only terms which tend to infinity as 
k z — 0. The amplitude coefficients of equations (58) also tend to infinity as k z — 0, with 
the result that every term in the stored energy expressions tends separately to infinity. 

In order to apply Rhodes' method to the cavity, one must distinguish between the l/k| 
type infinities and the infinities due to the amplitude coefficients. The reason for the 
infinite coefficients is that k z = 0 (or k z = mr/d) corresponds to the resonant frequency 
of the lossless cavity. A physically realizable cavity would necessarily have some ohmic 
losses so that these coefficients would not in fact become infinite. Alternatively, if the 
cavity were truly lossless, the aperture fields would have to vanish as k z — mr/d. On 
this basis, then, the infinities for which cancellation of terms are sought after the manner 
of Rhodes will be those infinities not caused by the amplitude coefficients. 

An examination of the expression for complex power at the aperture shows that it is 
infinite since the amplitude coefficients A mn and C mn become infinite as k z — 0. 

(See discussion following eqs. (62).) However, the infinities noted in the magnetic and 
electric volume integrals are not due entirely to the infinities caused by A mn and C mn , 
as discussed earlier in this section; they are also caused by the k| terms in the denom- 
inators resulting from the E z , H x , and Hy components. It is these latter contributions 
to the volume integral which must cancel identically when the volume integrals are differ- 
enced. Since the infinite complex power at the aperture is caused by the coefficients 
becoming infinite, the net infinite stored energy must result from these coefficients only'. 
Therefore, expressing the contributions which are infinite because of the l/k§ and the 
amplitude coefficients in terms of the contributions which depend on the amplitudes only 
will enable one to define "physically observable" stored energies for a shorted rectangular 
waveguide. 

By subtracting and adding the finite contributions in each volume integral and then 

grouping the negative finite contributions with the infinite contributions, each volume inte- 

* 

gral is then written as a new infinite contribution plus twice the finite contribution. This 
is equivalent to expressing the infinite contribution as a new infinite contribution plus the 
finite contribution. Thus, by following this procedure the two volume integrals are written 
as 
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Taking the difference in the two volume integrals (eqs. (58) and (60)) and regrouping terms 
results in the following: 



where ^i}, ^n}, {jni}, and {iv} are defined in appendix D. The terms {}} - k^{ll} 

and {ill} - {jv} are shown to vanish identically in appendix D. Thus, the magnetic 
and electric stored energies remaining are interpreted as "physically observable" by the 
following energy equations: 



III 


and 


«W e » v = ^ | 


< I I 2 fn e m| A mn | 2 + r m «n| c . 
m =0 n =0 
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(62b) 


It should be noted that these "physically observable" stored energies become infinite as 

k z — 0 |as well as for k z =~ where p = 1, 2, 3, . . ^ since the amplitude coefficients 

become infinite (see discussion following eqs. (57)). 

These "physically observable" stored energies can also be determined from the 
surface integral of the complex Poynting theorem as follows: 
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Dividing equation (65b) by 2oo yields the same expressions derived earlier for "physi- 
cally observable" stored energies. (The nonpropagating contributions are easily identi- 
fied with magnetic and electric stored energies since they are always positive. However, 
this is not true of the propagating contributions; identification of these energies is aided 
by examining which field components contribute.) 

Since the expressions given by equations (61) are to be used in determining quality 
factor Q, it is necessary that each be positive (or zero) at resonance. There is no prob- 
lem with the contribution from the nonpropagating waves since it is always positive. How- 
ever, the contribution from the propagating waves can be negative. But what must be 
considered, in general, is the case when both propagating and nonpropagating wave con- 
tributions exist simultaneously; the propagating waves can be limited by the choice of 
waveguide size, but the nonpropagating waves exist whenever a discontinuity is present 
in the waveguide, which is the usual case. Since discontinuities produce nonpropagating 
waves, the existence of propagating waves only can not occur. 

For the rectangular cavity-backed slot antenna considered in this paper, losses are 
inherently included due to the radiating slot. The net stored energy at resonance was 
always positive when a single propagating wave together with nonpropagating waves were 
assumed. This result was found to be true for many different size cavities. 

AMPLITUDE COEFFICIENTS FOR A NARROW SLOT BACKED 
BY A RECTANGULAR CAVITY IN WHICH A SINGLE 
PROPAGATING WAVE IS ASSUMED 

In this section the amplitude coefficients for a rectangular cavity which backs a 
narrow slot whose electric-field distribution is sinusoidal are determined. In order to 
determine the amplitude coefficients A mn and C m n for a narrow slot, the boundary 
condition at z = 6 (see fig. 1) must be applied. For an aperture opening of some kind 
at z = 6, the boundary conditions are applied to the transverse electric field. In the 
aperture it is assumed that the electric-field distributions are given by E xa and Eya. 
Therefore, at z = 6 (actually the boundary conditions are applied at z = 0 to simplify 
the derivation since later on 6 — 0 anyway), 


CO oo 

E xa= J I (A mn + B mn )cos^xsin^y (n * 0) (66a) 

m=0 n=0 

OO OO 

Eya = X X ( < " mn + ^mn)sin x cos y (m * 0) (66b) 

m=0 n=0 5 
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Multiplying by cos x sin y and sin x cos ^ y, respectively, and integrating 

a b 3- b 


over 0 to a and 0 to b, the coefficients become 
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(67a) 


(67b) 


(67c) 


(67d) 


Using the boundary condition at z = d for a shorted waveguide results in 
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where 



Therefore, equations (67) become 
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(68a) 
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(68b) 
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At the plane lim (z = 6), a narrow slot is located as shown in figure 1. A reason- 
6— 0 + 

able assumption for the electric-field distribution of such a slot is (for relatively deep 
cavities, ref. 7): 


Exa(y’.« = ^7 si" kff - |y 

E ya (y',k) = 0 




(69) 


Using the distributions from equations (69), the amplitude coefficients become 
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where 
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(propagating waves) 


(nonpropagating waves) 


Equations (70a) and (70c) can be integrated to give (see appendix E) 
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(71a) 


(71b) 
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Assume that only the m = 0 and n = 1 mode propagates, all others being evanescent 
(i.e., - < k < -V so that 

V b a / 
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n = 1, 3, 5, . . .J 
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Substituting the coefficients of equations (73) into the "physically observable" stored 
energy expressions given by energy equations (62) gives for «W m » v and «W e » v 
the following: 
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and 
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t=l x> 



Now, making use of the identities 


sin 2x = 2 cos x sin x 


and 


sinh 2x = 2 cosh x sinh x 
gives finally 
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( 77 ) 

(Equation continued on next page) 



In conclusion, equations (76) and (77) represent "physically observable" magnetic 
and electric stored energies for a section of rectangular waveguide which is shorted at 

depth d and bounded by a slot at lim (z = 6) in which the electric -field distribution 

6 - 0 + 

is assumed to be 



where <L is the length of the slot and a' is its width. The transverse dimensions of the 
rectangular waveguide are chosen such that only one field propagates (m = 0 and n = 1); 
all other fields are evanescent. 

RESULTS FOR THE RECTANGULAR CAVITY-BACKED SLOT ANTENNA 
IN WHICH A SINGLE PROPAGATING WAVE IS ASSUMED 

General 

When m = 0 and n = 1 are assumed for the propagating wave with all others 
taken as nonpropagating waves, the differences in the stored energies in the exterior 
region (for a slot, ref. 17) and in the interior region (eqs. (76) and (77)) are given, respec- 
tively, as 


43 









(Equation continued on next page) 
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(81) 


The total stored electric energy at resonance is needed for computing Q; i.e., the sum of 
electric stored energies from the half-space and the cavity regions. The half-space 
electric stored energy is given by Rhodes in reference 17 for a slot, and the cavity electric 
stored energy is given by equation (77). Therefore, 
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(Equation continued on next page) 
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For relative bandwidth the angular frequency derivative of equation (81) is needed, that is, 






(83) 

(Equation continued on next page) 
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The radiated power is given by (see ref. 17) 
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(84) 


The quality factor, relative bandwidth, and the input admittance given by equations (25), 
(21a), and (10), respectively, are repeated here for completeness 
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(85b) 



Approximate Solutions for Moderate Cavity Depths and Narrow Slots 

Since the summation on m in equations (81) and (83) is slowly convergent for 
a' « a, the numerical evaluation of the double summation in these equations is cumber- 
some. Adding further to the complications of evaluating these equations is the presence 
of the hyperbolic terms. The evaluation of equations (81) and (83) is greatly facilitated 
by making the following approximations: 

coth d 


coth d 

These approximations (eqs. (86)) place additional restrictions on the cavity size. From 
tables given in reference 20, the hyperbolic cotangent becomes very large for small argu- 
ment and approaches unity as the argument increases (around 3 or 4). Therefore, the 
approximation becomes less valid for decreasing argument. Physically, this means the 
deeper the cavity the more valid the approximation (coth ~ 1). This is also true for the 
assumed aperture distribution; that is, the sinusoidal distribution assumed in the slot is 
valid so long as the back wall of the cavity is not too close. 

By indicating a criterion that the hyperbolic cotangent functions must satisfy, the 
depth of the cavity is restricted, e.g., the arguments of equations (86) may be chosen so 
that 

d 


d 
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or 



(88a) 



2, 4, 6, . . A 

or 

1, 3,5,.. .) 


(88b) 


where in equation (88a) only n = 3 is necessary since for all other values of n the 
depth d is also greater. However, in equation (88b) the size of height a and width b 
will enter into determining which m and n must be used for d to be greater for all 
m and n. In addition to this restriction, the restriction resulting from the assumption 
of a single propagating wave must be satisfied; that is, 





or 
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J 

Even with the approximations (87) to (89) (coth ~ 1), the summation on m is still slowly 
converging for a' « a; hence, more approximations are needed to make the summation 
converge faster. (The details of the a* « a approximation are given in appendix F.) 
Finally, substituting the approximations into equations (81) to (83), yields the following 
solutions for moderate cavity depths and narrow slots: 
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(Equation continued on next page) 


51 



(Equation continued on next page) 


52 



53 



8|V r 


(2tt) Z 0 


£ k£< 
2_ 4 


Si(k£) - i Si(2k£) 
2 


sin k£ 


3 / 2 kt , 


Cin(k«) - - Cin(2k«) - ln(- 

2 \ 2kd. 


COS k£ + (92) 

M 


Input Admittance Calculations 

The input admittance of the rectangular cavity-backed slot is represented by 
Y = G + jB where the power loss in the cavity is assumed negligible. Explicitly, the 
conductance and the susceptance are given as 
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respectively, where equation (84) has been substituted into equation (85a). 

The input admittance can be represented as a function of any of the cavity or slot 
dimensions in wavelengths, as well as a function of frequency. For instance, a family of 
curves (G and B versus i/\) can be generated for fixed cavity size and varying slot 
widths in wavelengths. The representation of these curves is similar to the input imped- 
ance curves of the planar dipole antenna considered by Rhodes in reference 21 and of the 
monopole antenna considered by Jordan and Balmain in reference 5. Such curves would 
enable the user to design rectangular cavity-backed slot antennas. Once the operating 
frequency (resonant frequency, B = 0) is chosen, the physical slot length and width and 
the physical dimensions of the cavity can be determined from these curves. However, in 
order to compare such curves with measurements, either the frequency must be held con- 
stant and the physical length of slot allowed to vary while the remaining dimensions are 
held fixed or all dimensions except the slot length must vary in such a manner that the 
electrical lengths remain constant as the frequency is changed. Obviously, the first 
method of performing the measurements would be more desirable. 

Before presenting any design admittance curves, the input admittance represented 
by equations (93) and (94) is compared with experimental results (Long, ref. 22). In ref- 
erence 22, measured input impedance data for a slot (1 by 25 cm) backed by a rectangular 
cavity cross section (10 by 35 cm) with variable cavity depths over a frequency range of 
approximately 500 to 750 MHz are presented. Actually, only half the impedance was mea- 
sured since an imaging plane was used to bisect the slot lengthwise; the measured imped- 
ance must then be doubled. 

Admittance as a function of frequency using the representation given by equa- 
tions (93) and (94) for the same cavities and slot size with cavity depths of 13.395, 17.86, 
22.325, and 35.72 cm is presented in figure 3, and measured data are presented for com- 
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Admittance, millimhos Admittance, millimhos 


IIIMIIIHIIIIRII 


parison. The measured susceptance is in good agreement with calculations, particularly 
around resonance. The slope of the susceptance (which is related to bandwidth) at reso- 
nance becomes steeper as the depth of the cavity increases; the bandwidth becomes nar- 
rower. Also, the resonance frequency decreases as the depth of the cavity increases. 

The measured conductance is always greater than the calculated conductance; this may be 
caused by the losses in the cavity and ground plane, which have been neglected in the 
calculations. 


g ) Calculated 

• o ° G 1 

K f Measured 



Calculated 

Measured 



3 ) Calculated 





G I 

__ B Calculated 
° ° ® g | Measured 


3 ;?. 72 cm . 


1 J 


600 


Figure 3.- Input admittance as a function of frequency for a 1 = 1 cm, % = 25 cm, 
a = 10 cm, and b = 35 cm with varying cavity depths. 


56 



Since the representation of the input admittance has been shown to agree reasonably 
well with measured data, theoretical curves, based on this representation, were evaluated 
numerically and plotted as a function of slot length in wavelengths. In figure 4 the input 
admittance as a function of electrical slot length for three electrical slot widths (0.1, 0.01, 
and 0.001) are shown for increasing electrical cavity depths, respectively; the electrical 
cross section of cavity for these figures is rectangular (0.3 by 0.6). Similar curves for a 
square cavity cross section (0.6 by 0.6 wavelengths) are given in figure 5. Obviously, 
numerous sets of admittance curves for many different choices of slots and cavities could 
be computed. However, the calculated curves presented in figures 4 and 5 represent 
typical sets of curves that can be generated by equations (93) and (94). Such curves do 
offer invaluable design information; that is, once the operating frequency (resonant fre- 
quency, B = 0) is selected, the slot length and width and the cavity dimensions can be 
determined. Admittance curves as a function of frequency can also be generated for 
known physical dimensions, as was shown in figure 3. 



Figure 4.- Input admittance as a function of ) l/~K 
for a/A = 0.3 and. b/A = 0.6. 
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Figure 5*- Concluded. 


An examination of equation (94), which represents the susceptance as seen at the 
input terminals of the slot, shows that the contribution from the propagating wave (first 


term in eq. (94)) to this susceptance is zero whenever kd\l - — ) = — for odd values 


kb 2 


of n: the contribution of this term is infinite whenever kd 1 - / — = rur for all values 

kb/ 


of n. Its contribution to the total susceptance for other values of kd\ 1 - is finite 


and may be either inductive or capacitive; its contribution is inductive whenever 


nn < kd\ 1 - < (2n + 1) ^ for n = 0, 1, 2, . . . and is capacitive whenever 


(2n - 1) ^ < kd\ 1 - (— j < mr for n = 1, 2, 3, . . . . The contribution that the second 


term in equation (94) makes to the total susceptance is capacitive (this is the case for 
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n = 1 and m * 0 in the original double summation given by eq. (FI)). The contribution 
from the third term in equation (94) is inductive (this is the case for n > 1 in the sum- 
mations given by eq. (FI)), and the contribution from the external region, given by the 
last term, is either inductive or capacitive. 

As was noted earlier, the calculated conductance of the input admittance will prob- 
ably be slightly lower than its actual value since the analysis does not include cavity and 
ground plane losses. The calculated conductance, which is given by equation (93), is one- 
half the conductance of the slot when it is located in an infinite perfectly conducting ground 
plane and is free to radiate on both sides. 

Quality Factor and Relative Bandwidth Calculations 

Before the quality factor Q and relative bandwidth B.W. can be computed, the 
resonant frequency (or slot length) must be determined for a given set of cavity and slot 
dimensions (or slot width and cavity dimensions in wavelengths). The resonant slot length 
rather than the resonant frequency will be determined because it is more convenient to 
work with electrical lengths and because it gives a more general description of the slot 
length. This resonant slot length £ 0 /X i s found numerically by setting equation (90) to 
zero for a given set of electrical slot width a'/X and cavity size a/X, b/X, and d/X. 
The computer subroutine for determining the resonant lengths was verified graphically 
for the cross-over (B = 0) frequencies 620 MHz, 593 MHz, 574 MHz, and 530 MHz shown 
in figures 3(a), 3(b), 3(c), and 3(d), respectively. 

Once f Q /X is determined, its value along with the corresponding values of a'A, 
aA, b/X, and dA are used in the following equations to determine Q and B.W.: 


Q = 



(95) 


2P V 


B.W. = 


-u. 


-if 2w (. 

do>[ \ 


<< w m >> - «W e » 




u>=a> T 


(96) 
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where 



Si(k£) - ± Si(2k£) - — Cin(M) - Cin(2ki) - kH In 2 > cos k£ 
2 2 [ 
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(98) 

(Equation continued on next page) 



( 98 ) 

(Equation continued on next page) 


63 



- | Si(2k£) 


sin k£ + 


Cin(k£) - - Cin(2k£) 
2 


- In 



cos k£ + 



(98) 


and 


Pr = 


1 liJ C 


(2 if) 


2 Z 0 2 2 ^ 


Cin(kP) + 


Cin(k£) - ^ Cin(2k£) 

u 


cos k£ - 


Si(k£) - - Si(2k£) 
2 


sin k£ , 


(99) 

The quality factor Q and the reciprocal of the relative bandwidth l/B.W. using 
equations (103) and (104) were computed as a function of a'/X for a number of different 
size cavities. A typical set of such curves is shown in figure 6 in graphical form. 


From the curves shown in figure 6, it can be seen that both curves vary linearly 
with the log of a'/X; furthermore, their slopes are approximately the same. These char- 
acteristics were found to be true for the cavity sizes considered in this paper. Therefore, 


Q ~h log — + C! 




— - — ~ h log — + Co 
B.W. X 1 




( 100 ) 
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where h is the slope (which appears to be independent of cavity dimensions, h * -6) 
and Cj and C 2 are constants which depend on the cavity and slot size. These con- 
stants, which are the intercepts of Q and 1/B.W. axis, are quite large. However, for 
the cases considered in this research, the difference in 1/B.W. and Q is finite; that 
is, 


1 


B.W. 


- Q ~ C 2 - Cj — C 


( 101 ) 


where C is a finite constant which depends on the cavity and slot dimensions. There 
fore, 


1 

B.W. 


* Q + C 


( 102 ) 


The magnitude of 1/B.W. and Q determines the importance of the constant C; that 
is, if C is of the same order of magnitude as Q, then the reciprocal of relative band- 
width 1/B.W. is related to Q by equation (102). However, as Q becomes large, its 
relative magnitude compared with that of C is many orders of magnitude greater (C 
compared with Q is very small). Hence, 

11m rr-^Q (103) 

Q_»°c B.W. 

which is the asymptotic relationship usually assumed in the literature. It is concluded 
that according to the application of Rhodes' method of identifying "physically observable" 
stored energies the reciprocal relationship, usually associated with nonradiating systems, 
is equally valid for the rectangular cavity-backed slot antenna. 

Since Q and 1/B.W. were found to approximately satisfy equation (100) for the 
many cavity dimensions investigated, results are given in tabular form. Table I shows 
the computed resonant slot lengths for given electrical slot widths and cavity dimensions. 
The intercepts Cj and C 2 which are obtained from equation (100) are also given. 

From this table, therefore, the Q and the 1/B.W. (and hence, relative bandwidth) can 
be computed with very little effort. 
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TABLE I.- COMPUTED RESONANT SLOT LENGTHS FOR GIVEN 
ELECTRICAL SLOT WIDTHS AND DIFFERENT 
CAVITY DIMENSIONS 


1 


Q= -6 logi- + C i; 


1TW7 *_6 ] 0g aL + c 2 


r » 


£ 0 /X for a'/X 

of - 

C 1 

c 2 

10“ 2 10“ 3 lO -4 10' 5 

10 -6 

a/X = 0 . 3 ; b/X = 0.6 

- 4.5 

- 1.5 

0.500 0.500 0.500 0.500 

0.500 

- 3.2 

0 

.489 .493 .495 .495 

.495 

- 1.1 

2.0 

.479 .487 .496 .492 

.493 

1.0 

5.0 

.469 .480 .485 .488 

.490 


a/X = 0. 


a/X = 0. 


a/X = 0.6; b/X = 


a/X = 0.3; b/X 


0.497 

.483 

.467 


0.499 

.494 

.489 



- 3.5 

- 1.0 

0.493 

0.495 

0.496 

0.497 

0.498 

- 2.7 

-.3 

.484 

.490 

.493 

.494 

.495 

- 1.5 

1.0 

.478 

.486 

.489 

.491 

.493 

.5 

3.5 

.470 

.481 

.486 

.488 

.490 


- 2.7 

- 0.7 

0.489 

0.493 

0.494 

0.495 

0.495 

- 2.0 

.5 

.482 

.488 

.491 

.493 

.494 

- 1.0 

1.5 

.477 

.485 

.489 

.490 

.492 

.5 

3.0 

.471 

.481 

.486 

.488 

.490 


- 2.7 

- 0.5 

0.487 

0.492 

0.494 

0.496 

0.497 

- 2.0 

.3 

.480 

.481 

.490 

.492 

.494 

- 1.5 

.8 

.475 

.484 

.488 

.491 

.492 

-.4 

2.5 

.471 

.481 

.486 

.488 

.491 


0.500 

.497 

.494 

.491 

0.498 

.495 

.494 

.491 


0.497 

.494 

.493 

.492 


0.497 

.494 

.493 

.492 



0.500 

.497 

.495 

.493 


0.498 

.495 

.494 

.493 


0.497 

.495 

.494 

.493 


0.498 

.495 

.494 

.493 


0.499 

.495 

.492 




TABLE I.- COMPUTED RESONANT SLOT LENGTHS FOR GIVEN 
ELECTRICAL SLOT WIDTHS AND DIFFERENT 
CAVITY DIMENSIONS - Concluded 

Q = -6 1og^ + C i; . * -6 log ^- + C 2 


d/X 

C 1 

c 2 



W x 

for a'/X 

of - 



10 -2 

1— i 

o 

i 

00 

10-4 

o 

1 

CJI 

10“ 6 

10 -7 

10' 8 





a/X = 0.3; b/X = 

0.8 




0.3 

-2.0 

-0.4 

0.496 

0.497 

0.498 

0.498 

0.498 

0.498 

0.498 


-.5 

2.0 

.419 

.486 

.489 

.491 

.493 

.493 

.495 

.5 

1.5 

8.0 

.457 

.471 

.478 

.482 

.485 

.487 

.489 





a/X = 0.4; b/X = 

0.7 




0.3 

-2.0 

0.0 

0.491 

0.494 

0.495 

0.496 

0.497 

0.497 

0.497 

.4 

-1.0 

1.2 

.480 

.487 

.491 

.492 

.493 

.494 

.495 

.5 

.0 

3.3 

.469 

.480 

.485 

.488 

.490 

.491 

.492 





a/X = 0.4; b/X = 

0.8 




0.3 

! -1.5 

-0.5 

1 0.490 

0.493 

0.494 

0.496 

0.496 

0.497 

0.497 

.4 

-1.0 

1.2 

.478 

.489 

.489 

.491 

.492 

.493 

.494 

.5 

j 1.0 

5.5 

.461 

.474 

.481 

.485 

.487 

.489 

.490 


CONCLUSIONS 

It is concluded that the time-average "physically observable" stored energies in an 
infinite rectangular waveguide and in a rectangular cavity can be expressed in a form 
which requires only the transverse electric and longitudinal magnetic field components. 
This same conclusion had been reached earlier for the planar aperture antenna. For the 
infinite waveguide case, the propagating fields do not contribute to these stored energies 
whereas the nonpropagating fields do. However, both propagating and nonpropagating 
fields were found to contribute to the stored energies for the rectangular cavity case. 

The representations of these energies were found to be analogous to the representations 
given for the planar aperture antenna; although, in the cavity case the representations are 
modified by the presence of the propagating field contributions and functions which depend 
on the cavity depths. It is concluded that the method of Rhodes can be adapted to the 
infinite waveguide and finite waveguide (cavity), etc. 
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Input susceptance calculations were found to agree with measurements, particularly 
around resonance. The slope of the susceptance at resonance became steeper as the 
depth of the cavity increased; the bandwidth became narrower. Also, the resonant fre- 
quency decreased as the depth of the cavity increased. The calculated input conductance 
was always less than the measured input conductance. This difference was attributed 
to neglecting the losses in the walls of the cavity and ground plane in the analysis. 

The asymptotic relationship which is usually assumed to exist between the quality 
factor Q and the reciprocal of relative bandwidth for nonradiating systems was found 
to be equally applicable for the rectangular cavity-backed slot antenna. For lower values 
of Q, the reciprocal of relative bandwidth and Q are approximately related by constants 
which depend on the slot and cavity dimensions. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 31, 1975. 
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APPENDIX A 


DERIVATION OF FIELDS INSIDE A RECTANGULAR WAVEGUIDE SECTION 

From Maxwell's equations, assuming an ej wt time convention, a wave equation is 
obtained in terms of the electric field E, 

V 2 E + k 2 E = 0 (Al) 


where k = Expressing E in its rectangular components, equation (Al) for the 

E x and Ey components becomes 


9 2 

— E x (x,y,z,k) 
9x 2 


2 2 

+ E x (x,y,z,k) + - — E x (x,y,z,k) + k 2 E x (x,y,z,k) = 0 

9y“ 9z^ 


(A2a) 


2 2 q2 

~2 E y (x,y,z,k) + -£-5 E y (x,y,z,k) + E y (x,y,z,k) + k 2 E y (x,y,z,k) = 0 (A2b) 

9x^ 9y 9z 


By separation of variables, a solution to equations (A2) can be written as 


E x (x,y,z,k) = (A cos k x x + B sin k x x)(C cos kyy + D sin kyy)^Fe ^ Z + Ge^ zZ j (A3a) 

Ey(x,y,z,k) = (A' cos k x x + B' sin k x x)(C' cos kyy + D' sin kyy)^F'e ^ + G'e^ j (A3b) 
2 2 2 9 

where k x + ky + k z = k^ and A, B, etc., and A', B', etc., are arbitrary constants. 

The coefficients will now be determined for a rectangular waveguide (see fig. 1). 
The boundary conditions on Ey at x = 0,a are 

E y = 0 (A4) 

that is, for all y and z within the cavity, 


A'(C' cos k y y + D' sin kyy^F'e "'^ 22 + G'e jkzZ ) = 0 (A5a) 

(A' cos k x a + B' sin k x a)(C' cos k y y + D' sin kyy)^F'e ^ + G'e^ ) = 0 (A5b) 
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APPENDIX A — Continued 


Equation (A5a) implies that A' = 0. Hence from equation (A5b) for B' * 0, it follows 
that sin k x a = 0, or that 

k x a = mir (m = 0, 1, 2, . . .) (A6) 

Therefore, equation (A3b) becomes 

E y = B' sin x(C* cos k y y + D’ sin kyy^F'e"^ 2 + G'e^ 2 ) (A7) 

Applying the boundary condition E x = 0 at y = 0,b yields 

E x = (A cos k x x + B sin k x x)D sin — y^Fe ^ zZ + Ge^ zZ j (A8) 

where ky was found to have the value nn/b for n = 0, 1, 2, . . . . Rewriting equa- 
tions (A7) and (A8) with k x = and ky = 


E x = d(a cos HI x + B sin HI x)sin 'f y(Fe' iM + Ge )kzZ ) (A9a) 

Ey = B’ sin HI xfc’ cos ~ y + D' sin ~ yj^F'e~^ kzZ + G'e jkzZ j (A9b) 


At this point Maxwell's divergence equation V • D = 0 is used to determine the remain- 
ing electric field component E z 


V • D = 0 


3E X 9Ey 

+ + 

3x 9y 


9E Z 

9z 


= 0 


(A10) 

(All) 


9E Z 

9z 


9E X 9E^ 


3x 


8 y 


-D sin M x + B cos ^ I x)sin f y^Fe" jRzZ + Ge jkz j 


. B-HsinMx 
b 


“ jkzZ + Ge jkzZ ) 

C* sin ^ y + D’ cos ~ yj^F’e ^ kzZ + G’e^ kzZ ) (A12) 
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APPENDIX A - Continued 


Integrating on z, equation (A12) becomes 


E z = -D HZ [_A sin HI x + B cos ^ xjsin ^ y|-±— e ' + 7 £- e 


HI xjsin HI y f-2L e‘ jkzZ + JkzZ ' 

a / b l-jkz jk 2 


. B ’ H sin HI x(-C' sin HI y + D’ cos HI y )(Zl e'^ 2 + e* 22 ' 

b a l b b y hk z jkz 


(A13) 


The constant of integration has been set equal to zero since it would not represent a 
traveling wave in the z direction. Applying the boundary conditions that E z = 0 for 
x = 0,a and y = 0,b, 


■D SSI B sin 2S.I -X- e' ikzZ + e illzZ ) + 0 = 0 

a b l-]k z jk z / 


(x = 0) 


which implies B = 0 and 


0 - B’D' 22 sin — xf-~ e‘ ifez + e j H = 0 
b a V - jk z jk z / 


(y = o) 


(A14) 


(A 15) 


which implies D f = 0. The boundary conditions at x = a and y = b give redundant 
information. Substituting these conditions into the equation for the electric field compo- 
nents yields 


E x = DA cos 222 sin 21 y( F e' itzZ + Ge iM 

a fc) 


E v = B’C’ sin HH cos — ylF'e"^ 2 + G'e'’ kzZ 
b b 


E Z = DA HI sin HI x sin HI y / JL e ' jkz2 + e jkzZ ' 

a a b J (jk z jk z 


+ B'C HI S i n HI x sin HI y [ JLL e " jkzZ + JS1 e jkzZ 

b a b \-]k z jk z 


(A 16a) 
(A 16b) 


(A16c) 
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APPENDIX A - Continued 


Simplifying equations (A16) with the substitutions DAF = A m n, DAG = B mn , 

B’C’F* = C mn , and B'C'G' = D mn and recognizing the dependence of E x , Ey, and E z 
upon m and n these quantities (eqs. (A16)) may be written as 


fe) mn " [ A mn e ■* Z + 


B mn e jkz jcos^xsinM 


(Al7a) 


(®y) mn " fcmn e + D mn e^ z J sin x cos ^ y 


(A17b) 


(E z ) 


mn 


"i k z 


a + r )p”^ kzZ | m7r R 4. n7r n \J kzZ 

a A mn + — L mn le - I— Bmn + -g- D mn le 


fm77 a , nir 


sin x sin ^ y 
a b J 


(A17c) 


Thus the most general solution of equation (A3) will be given by 


OO oO 

E x = Yj /L ® x ^mn (A 18a) 

m=0 n=0 

OO OO 

E y= ^ ^ ( E y)mn (A18b) 

m=0 n=0 

OO OO 

E 2 = l X ^ Ez * mn (A18c) 

m=0 n=0 

From Maxwell's curl equation 

V X E = -jcojiiH (A 19) 

the magnetic field can be determined. Hence, the magnetic field components are deter- 
mined as follows: 
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( H y)mn " 


APPENDIX A - Continued 


1 

9(E Z ) 
v z, mn 

8 <%>nm 

ja)JU 

9y 

3z 


1 Jmr/b 


jw/x Vjkz 


a . 321 c i p 
— A mn + 1“ L mn e 


■jkz z / m7r rj n-ff p. \ Jk z z 

“ I “mn + ■jj" u mnl e 


x sin ™ x cos ^ y - f-jkgCmn e ^ zZ + jk z D mn sin — x cos ~ y^ 


— — 

jw/i\-jk Z/ 


m 7 T\/n 7 r 


a b 


Amn + 


fr*4 


Cmn'* e 


-jl^z 


mir\nir\-o 
~ ] b" ,Bmn 




D mn^ e 


jk z z 


sin x cos — y 
a b 


k z w/x 


— F )A mn + 


H\ + k 2 

b 1 +kz 


^mn / e 


-jk z z 


rmr\fniT\jy 

r b- ,Bmn 


?) +k i 


^mn/ e 


jk z z 


sin™xcos^ y 

a b 


(A20) 


i 

' 8 < E *>mn 

«Wnm 

jwp 

| 8 z 

dx 


1 


jk z A mn e“ jkzZ + jk z B mn e jkzZ ) cos S 2 x sin 22 y - (- 1 


jkz i 


mi t 


Mnn 


(A 21 ) 

(Equation continued on next page) 
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APPENDIX A - Continued 


f C mn je' j 

/fr 


' jkzZ . (m Bmn + f) Dmn e jkzZ 


“IcosMxsin^y; 
a a b ' 




k 2 + S£ 


A m „ * |!Sa(ffi!)c mn >e' ilIzZ 


*♦(?) 


B 


mn 


. ImnXfnn^ V j k z z nur „ . nn „ 
‘'TtN 1 cos— xsm-y 


(Hz) 


mn 


1 

9 < E y>mn 


9x 

1 

safe 


~rmn 
a \ 


'mn 


a y 


-jk z z Jk z z 


2221 x cos — y - — ( 
a b b \ 


-jk z z 


jk z z] mir nir „ 

+ Bnin e ) cos x cos — v 

/ b a 


1 

jW|U 


'mn n nn * V“jk z z 
L mn - -g- A mnl e 


j- / mir n - 077 R U k z z 

+ 1 -^— u mn ■ -g - B mnl e 


x cos 2221 x cos 221 y 
a b J 


(A21) 


(A 22 ) 
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APPENDIX A - Concluded 


Hence, 


oo w 

H X = I ^ Wmn 


m=0 n=0 


oo oo 

H y = Z Z ( H y)mn 


m=0 n=0 


H Z -I ^ (H z ) mn 

m=0 n=0 


(A23a) 


(A23b) 


(A23c) 



APPENDIX B 


PROOF THAT THE CONTRIBUTIONS FROM THE PROPAGATING 
WAVES TO THE STORED ENERGIES CANCEL WHEN 
THE VOLUME INTEGRALS ARE DIFFERENCED 

That the contributions from the propagating waves to the stored energies cancel 
when the volume integrals are differenced is proved as follows: 

From equations (45) and (46), 


lim 


00 00 

r 


r -1 


2 

r -1 

ab^x y y 1 j 

^m e n 
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\b J ^ 
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16 ^ 0 2 2] 

m=0 n=0 w M 
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K z 
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L mn -j^“ A mn 
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d -°° 16 co 2 ju 2 
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- 2k 
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APPENDIX B — Continued 


where 


4 - * 2 - (¥) 


'nir\ 2 

, b J 


- k 2 
z2 K 


k 2 = k 2 [—} (—) 

z ~ [a J ~ [b J 


J 


Substituting equations (B2) into the right-hand side of equation (Bl) yields 


(B2) 


ita*y y -i-„ 

d — 00 16 J 2 

m=0 n=0 w ^ 


r 




2 

r 

y € 
A m n 

MHa + 

(?f ♦ 4 

p 

y € 

d + n m 

k 2 /mn\ 2 

[4 

UA>>r mn 

^mn 

° 2 

Z l~j 


l mn 




"mn 


d + e m e n 


m7r ^ n?r A 
L mn " A mn 


00 00 


> - Um ^ y > 

I 16 A 


/ 


m=0 n=0 


Vn e 


n m 


■A-mnl d 
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Rewriting equation (B3) yields 
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Hence, 


lim 


00 00 



r 


2 

" " 

abii y y 1 J 

V € 
r m n 

/W\/W\a 

/sif + fe2 

r* 

d 1 n m 

k 2 + /sif 

16 L, Li 2 ., 2 1 

m =0 n =0 w M 


UA b r mn 

l b / ^ 

L mn 

+ ^ 

Z \ a / 


Mnn 


m ? t 


f n Ar 

\b J mn 


d + 

mn 


r _ lLn a 
— c mn -jj- A mn 


oo 00/ g 

d f " TJ ^ £ Kn e m[ A mn| d 


m =0 n =0 


+ y m e n 


' mn 


d + 



m 7 r 

k 2 

K z 

a 


A + — C 
A mn + ^ L mn 


dUo 


(B 5 ) 
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APPENDIX C 


PROOF THAT THE BRACED TERMS GIVEN IN EQUATION (49) CANCEL 


That the braced terms given in equation (49) cancel is proved as follows: 
From equation (49) 




(Equation continued on next page) 
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- lim 


abe V \ 

1 ^ *\j c 1 A 

2 

— •w C 

n 

2 y m y n 

mu 

16 A A 

m=0 n=0 

\ ^n fc m | A mn 

y m e n 

l mn 

k 2 

K z 

a 


a 4. n7T r 
A mn + "g" L mn 



(C5) 
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APPENDIX D 


PROOF THAT TERMS IN EQUATION (61) CANCEL 
In order to show that all terms in equation (61) cancel, first let 





M8 


APPENDIX D - Continued 





(D3) 

(Equation continued on next page) 
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Substituting equations (D6) into the right-hand side of equation (D5) yields 
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(D8) 

(Equation continued on next page) 
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(DIO) 

(Equation continued on next page) 
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■■Ill ' 



Then, from equation (61), 



Since the term 




in equation (D12) is identical to the 


same term in equation (D4), it follows that 
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APPENDIX E 


EVALUATION OF THE INTEGRALS SHOWN IN EQUATIONS (70) 


Consider the integral given in equation (70a) as 


f ,/2 f' /2 

d y'= -St/2 d x'= -a'/2 a l 2 


y’ 


sin dx’ dy' 
b J 


f*/2 / I 

= V n \ sin — cosky' 

oJ -f/2V 2 1 


IrjO 

cos — sin k 
2 


y 


s in ™( y . + b 

b r 2 


dy’ 


f*/2 / 

= V„ \ sin — cos k 

°J.£/2\ 2 


y’| - cos ^ sin k|y’|j^sin cos ^ + cos ^2^ sin £pjdy’ 

V Q C (sin ^ cos ky’ + cos — sin ky*j (sin cos — + cos M2L sin — Wy’ 

0j -£/2\ 2 2 V\ b 2 b 2 / 

+ f (sin ^ cos ky’ - cos ^ sin ky’Ysin cos ^ + CO s 2^2- sin — ' )dy’ (El) 

«^0 V 2 2 A b 2 b 2 / 


Now let y’ = -y’ in the first integral so that 


f * /2 r ’ /2 ^ sl „4 

T=-l/2 J x--a'/2 a ’ \2 


= V 0 f (sin ^ cos ky' - cos ^ sin ky’V-sin n7 ^ L cos — + cos M2L S j n EE 
*J0\2 2 yyb 2 b 2 y 


y’ 


sin dx' dy* 
b J 


+ (sin — cos ky' - cos — sin ky'Vsin EEL. cos — + cos S i n HE 
y 2 2 J\ b 2 b 2 , 


dy’ 


= V n f (-sin ^ cos ^ cos ky' sin — ^ + sin ^ sin ^ cos ky' cos HE2L 
u Joy22 b 22 b 

+ cos — cos — sin ky' sin - n77 ^ - cos — sin — sin ky' cos EE L. + s in — cos — (E2) 

22 b 22 b 22 


(Equation continued on next page) 
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APPENDIX E - Continued 

Making use of equations (E3) and (E4) combined with equation (E2) results in 


f V2 f' /2 ^sinkff-ly'LinSHdx’dy' 
Jy»=_£/2 ^x^-a'/Z a \ 2 ' '/ b 


V Q sinf 






sin — ^ 

M 

k 

1 


l L 


sinlk - 2ai + sinjk + H)i 

b 12 V b 12 


nir 


b/2 


sinfk - — - sinfk + — 


b/2 


+ cos ^<k| 


cos[k- f f + cosk tf | 


m t 


co S (k - - cos(k 


- 2k 


(E5) 


Using trigonometric identities, 


f 1/2 f aV2 ^sink/i- y-'lsin^dx'dy' 
^y'= -^/2 ^x'= -a’/2 a \ 2 / b 


V Q sin 2^ 
0 2 


V _ I _ 
k ' b 


kfsin 2 — cos — - + cos 2 — cos — -) 
\ 2 b 2 2 b 2/ 


— fsin — cos — sin — — 
b 2 2 b 2 


sin — cos — sin — - k cos — 

2 2 b 2 


V n sin — / \ 

— 2|k cosHL k cos M) 


T 2 /W 

k ‘ T 


b 2 


(E6) 


or for n even, 


nl/2 C a /2 X°sinkf|- y’l) sin 2221 dx*dy' = 0 
Jy'=-£/2 ^x’=-a'/2 a \ 2 '/ b 


(E7a) 
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or for n odd, 


r £/2 r ^ X 2 sin kfl - I y’f)sin ME dx’ dy’ = (-1) 2 2V 0 k 
«Wf/2 *Wa’/2 a’ 2 l y l b 0 


n_ ^ COS — — - COS — 

b 2 2 


k 2 - m 


J y'= -f/2 J x’=-a’/2 


b q 

Now consider the integral given in equation (70c) with y = — + y’ and x = — + x' 

Li Li 


(E7b) 


£/2 r a'/2 V Q /« | 

\ -£sink|- y’ 

J y’=-£/2 J x’=-a’/2 a V 2 


f 


cos sin M2, dx* dy’ 
a b 


= C £ /2 r /2 Z 2 sin kfl - I y’ )cos ^M/| + x^sin M/k + y* 1 ] dx» dy’ 

‘ V =-£/2 *' x *=- a ’/2 a \ 2 / a \ 2 / b \ 2 / 


V 0 II k£ 

— f- \ sin ^ cos k y f - cos sin k 

a f V=-*/2\ 2 11 2 

x r V2 /cos 22 cos SMHE! _ sin 222 sin dx’ dy ’ 

Jv’=-n»/ 5 ! I 2 a 2 a / 


' . nir n7ry' nir . n7ry*’ 

sin ^ cos — + cos — sin — J — , 
1 2 b 2 b 


Since 


r a’/2 

\ dx’ = 

J x’=-a’/2 


m77 . 

pos sin 

m77X' 

V/UD _ OIAi 

z 

a 

m7r 


a 



a’/2 


■a’/2 


. m7r 

n poc 

m7TX 

oill UJd 

u 

a 

m7T 


a 



a’/2 


■a’/2 


. r, mir mira’ 

r a /2 dx , _ 2 cos — sm jr 

*^x’= -a’/2 222 

a 


(E8) 


(E9) 
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APPENDIX E - Concluded 
Hence, equation (E8) can be rewritten as 


nSL/2 r>a'/2 v, 
-jf/2 ^x’= -a'/2 a 


° sin - j y' jjcos 21p£ sin 22X dx' dy' 


„ m 17 . m 7 ra' nir £ k£ 

2 cos — — sin — — v cos — - - cos — 

2 2 a v o ou mr b 2 2 

2 k sin — 

tmt a' 2 „ , ^2 

— k 2 - /EI 


(E10) 


f‘ /2 f’ /2 X2 slnk (| 

J y'=-£/2 J x'=-a'/2 a \ 2 


y’ 


) cos 21H sin 22Z dx' dy' = 0 
1 a b J 


(n even or m odd) 


(Ella) 


f i/2 f aV2 — sin kfl - 
J y'=-&/2 J x’=-a'/2 a* \2 


y' 


| cos — , — sin 2£Z dx* dy* 
1 a b J 


Sil rn sin I^Sl cos H | - cos M 

= (- 1) * (- 1) 4V 0 k (n odd or m even) 

a* S3JL . / ,2 

“ * 2 -(f) 

(El lb) 

where the y* integration has been evaluated as in equation (E8). 


98 



APPENDIX F 


APPROXIMATION FOR a’ « a IN EQUATIONS (81), (82), AND (83) 

Rewrite the summation terms given in the expression for the difference in stored 
energies beginning with equation (81) (coth ~ 1) as 


r 

n=3,5 


n7rJ ^ k£ 

Icos cos 

V 2b 2 

3/2 


+ 2 


MY . k 2 


00 oo 

1 l 




m- 2; 4 n=l,3 ^ Sg. j | |^f + . ,2 (f f - k* 

3 y l b j \ / 


oo 

- I 


(cos - cos M) 2 

\ 2b 2 ) 

1 

(mf . k 2 

ll/mf . k 2 


l b 


Li 



/if . k 2 rn=2\ ^ 


'ssf + /ir -k 2 


ay 'lb/ 


(FI) 


Write the function 


as 


fsin^ 
2a 

mra* 

2a 


in its integral form (the integral from whence it came) 


/sm saa’ 

2a 
mra’ 


= f cos £121 (| + x) dx C ^ cos -S^(§ + x ,N ) dx' 
a ’ J -a*/2 a \ 2 / a J -a’/2 a \ 2 / 


2a 


-i- f ^ r cos W§ + x) cos — (§ + x ,N ) dx dx’ 
( a if *^-a'/2 -a'/2 a \ 2 / a \ 2 / 


(F2) 
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Let x = and x’ = SS so that 
2 2 ’ 


/sin 


mira 


mna. 

2a 


^r~ ) = t f f cos ~ir~ ( a + a ' x ) cos ( a + a'x') dx dx’ 
* / 4 2a 2a 


mir , 


(F3) 


Expand the cosine product 


cos si (a + a’x)cos SI (a + a’x’) = cos SI CO s SZS^E . S in S2 sin SI^ cos SI 
2a ’ 2a \ 2 2a 2 2a /\ 2 


x cos SI*'*! - S i n BE sin m7ra ’ x ’' 


2a 


2a 


cos 2 HI cos Sn5E COS - sin HI CO s HI 

2 2a 2a 2 2 


X Us SHS5E sin + cos =- sin 

\ 2a 2a 2a 2a 


+ sin2SIsinSI a ^ si nSI a ^ 

2 2a 2a 


(F4) 


Substituting equation (F4) into equation (F3) gives 


sin SS’ 
2a 

mra.’ 

2a 


• if 1 f 1 cos2Hl 

4 J_ i J_i 2 


cos SlI'x cos ra7ra ' x 


2a 


2a 


dx dy’ 


• if 1 f 1 

4 J-l J-l 


_ m7ra f x mm'x , , f 

cos — cos dx dx T 


2a 


2a 


1 pi 


-n 


cos m7ra x cos m7ra x> dx dx’ 


0 ^0 


2a 


2a 


(m even) 


(F5) 
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APPENDIX F — Continued 

Consider the summation in which n = 1 (last term in equation (FI)) 




The summation on the right-hand side of equation (F6) is approximated by replacing 



The summation on the right-hand side of equation (F7) can be summed in closed form by 
rewriting as 


mira'x' 


°E, cos ( x + x') + cos n y ra< (x - x’) 
\ 1 23 . 23 

Lj 2 m 

m=2,4 




(F8) 
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APPENDIX F - Continued 


Consider the following sums (see ref. 19): 


Jmx 


m 


= -In 


m=l 


(l - e*) 


and 


* P jm(x+7T) 

e = -In 1 + e JX 


I 


m 


m=l 


where the last and first are written, respectively, as 


Jm(x+7r) 


OO OO 

jmx v 


m 


m=l 


I ^ + l 

m=l,3 m=2,4 


+ / = -In 1 + e^ x 

m Lj m 


and 


jmx 


OO 00.00 

y y e^f + y 

Zv m L m Li m 

m=l m=l,3 m=2,4 


= -In 1 - e* 


Now, add the two series of equations (F10) 


m=2,4 


jmx 


m 


= -In (1 + e^ x |fl - e'* x )= -In 


£ 


£ 


2e 2 cos ^ 2(-j)e 2 sin ^ 


In (cos - sin -] + 2 In 2 + j(x - ^ 

l 2 2 j \ 2j 


lnh| sin xj + 2 In 2 + j(x 


00 

1 


m=2,4 


e ]mx _ 1 

m " 2 


In sin x + ln2 + j(x-- 

2 


Y' p jmx . 

Re ) — = - — (In sin x + In 2) 

Li m 2 

m=2,4 


(F9a) 


(F9b) 


(FlOa) 


(FlOb) 



(Fll) 

(F12) 

(F13) 


102 



APPENDIX F - Continued 


Therefore, replacing x 


by 


i 1 


and 



- x’), equation (F8) becomes 


“ cos cos mraV 


V ” 2a 
m=2,4 


rmr 

a 


2a _ a J ^ 

27T V 2 


In sin (x + x’) + In 2 

LkC l 


In sin (x - x’) + In 2 



s- 

r - n 

- 

1! 

2 In 2 + In 

t 

sin •2^-(x + x f ) 
2a 

+ In 


sin -22- (x t x r ) 
2a 


_a 

47t| 


2 In 2 + In — (cos - cos 

2 \ a a , 


For small values of y Icos y ~ 1 - ), 


oo 


I 

m=2,4 


cos 


m7ra'x 

2a 


cos 


mTra'x' 

2a 


BE a'-0 

a 




Therefore, 


2 



00 

1 


sin 


mjnr 

2a 


m=2,4 


mwa' 

2a 



(F14) 


(FI 5) 


(F16) 

(Equation continued on next page) 
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a'— 0 



In 2 - 3 


Next, consider the case when n > 1; the summation on m is 


2 


OO 


l 

m=2,4 


'sin 


niTra* 

2a 


m7ra’ 

2a 


\ 2 

/ 



= 2 



m7ra’x m7ix f 

cos — cos 


2a 


2a 


dx dx f 


mi t 
a 


n7r 


- k 



OO 


l 


cos 


niTra* 

2a 


(x + X*) + cos 


m7ra* 

2a 


(x - x’) 



dx dx' 


The approximation made for the case when n = 1 is not valid here since (-2^ 


always greater than 


Mf _ k 2 


; hence, other approximations must be sought. 


(F16) 


(FI 7) 

2 

is not 
In order 


104 


APPENDIX F - Continued 


to convert this summation to a more rapidly converging series, use is made of the 
Fourier transform (ref. 19): 


Jo)(y+y’) r , 

e dw = 2 Kq |o;(y + y') 

^U> 2 + a 2 


(F18) 


where KQ(ay) is the modified Bessel function of the second kind, and of the definition of 
the Poisson summation formula 


i f(2ma> 

_ oo ' > m=- 00 


(F19) 


where f(y) is the Fourier transform of S(a>). (See ref. 19.) Convert the series to be 
summed by means of Euler's identity as 


OO 


I 

m=2,4 


cos n l 7r ~ (x + x) + cos 
2a 


m7ra’ 

2a 


(x - x’) 



OO 

\ l 

m=±2,±4 


.mra/ 

e 2a 


(x+x') 


+ e 


jssW) 




00 

1 

m=0,±2 


.jinn-a' 
e 2a 


(x+x') 


• imra' 
+ e 2a 


(x-x) 




(F20) 


Then from the Poisson summation formula 


I 



00 

IT X + 2ma ) 

- 00 


(F21) 
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APPENDIX F - Continued 
By replacing a with a/2, even terms are obtained thus ly 


,2m7ry 


I e 

m=0,±l 


■ 2 


.niry 


2rmr\ 2 + a 2 n=0 ’ ±2 


= ^ ^ K^ay + ma) (F22) 


(f) + 


Applying this result to the case at hand, equation (F20) becomes 



The modified Bessel function of the second kind decays rapidly as the argument increases 
therefore, m = 0 should be sufficient so that 



(F24) 
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APPENDIX F - Continued 


Hence, equation (F17) becomes 




For small argument, 



(F25) 




(F26) 

(Equation continued on next page) 
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(F26) 


where y is Euler's constant = 0.5772157. Therefore, the summation terms in the 
expressions for the difference in stored energies become 
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(F27) 

(Equation continued on next page) 



APPENDIX F - Continued 



Hence, equation (81) is now written as 


/ \ 4 v o 

2^(«W m » - «W e »J 


4 

v ol 

2 “l 

If i rt ktV 

l^cos 55 - cos T j 

abZ 0 | 

l 

fo / \ 2 ~ 
k 2 - fs\ 

\ h ) 

3/2 


cot 


k 2 - 


gfd 


TT G lr 0 

+ al C0S 2b- C0S 1T 

+ 7 r 

* 2 - a 


l„|Sl + ln2-| 


I 1 


oo nirti. k£ 

Icos 2b- - cos T 

+ - / -> - y - In 


n=3,5 (rvrr^ _ ^2 


F ' “ k2 


8 V f 

— <s MM) + 

(2tt) 2 Z 0 4 


Si(k£) - | &(2k£) 

a 


cos k£ 


3/2 

Cin(k£) - i Cin(2k£) - In - - - 
2 2a 


sin k£/ 


(F27) 



(F28) 
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Using the same approximation, the representation for 2co («W e ») 


(eq. (82)) is 



X <2 ' 7 ‘ ln 




^-k 2 d^ 
b j 2 


7 


01 f[Si(k£) -M C in(kf) 


(2tt) 2 Z 0 8 


(k£ - sin kH) ln ~ + jsi(k£) - | Si(2kf) - Mjcin(kf) - Cin(2kf)] 


- k£ In . 2^ cos k£ + jcin(k£) - ± Cin(2kf) + M jsi(kf) - Si(2k£)] 



Equation (83) is therefore written as 


(F29) 


<x> 


_d_ 

dw 



» - «W e ») 
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(F30) 

(Equation continued on next page) 



APPENDIX F - Continued 



(F30) 

(Equation continued on next page) 
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APPENDIX F - Continued 




Cin(k£) - - Cin(2k£) - In 
2 



x cos k£ + 



Rewrite equations (F28), (F29), and (F30) as 


(F30) 



1 e^/^k£ 

+ Cin(k£) - - Ciri(2kj?) - In - ^ 

c* ZK3, 


sin k£ 
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L K ,N 

+ In 2 - - 

\ka , 

/ 2 


ka V 
It Li 


oo / n?rk£ 

'cos — — — cos 


2kb 


¥) 


n=3,5 


ni ) 2 - 1 

kb 


x 


/f-r-in 


mr\ 1 ka| 
kb/ ~ 2 


A 


/ 


— f fsi(kjg) - M Cin(k£) 
(2 W ) 2 Z 0 8\ 2 


+ (k£ - sin k£) In + jsi(k£) - ~ Si(2k£) - M[cin(k£) - Cin(2kje)J 
- k£ In 2^ cos M + jcin(k£) - | Cin(2k£) + ^ j^Si(k£) - Si(2k£)] 


!} s “ 


+ In — > sin kjp) 


(F32) 



(F33) 

(Equation continued on next page) 
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(F33) 

(Equation continued on next page) 


APPENDIX F - Concluded 




- - Si(2k£) 
2 


sin k£ + 


Cin(k£) 


- | Cin(2k£) - In 



x cos k£ + 



(F33) 
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